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ABSTRACT 

Most analytic and computer simulation models of the Navy's 
CONREP phase of underway replenishment operations have assumed expo- 
nentially distributed service times for two reasons: (1) lack of 

better estimate of the true distribution, and (2) such an assumption 
leads to more tractable mathematical solutions. 

By analyzing several different combinations of replenishment 
vessels and combatant ships, it is demonstrated through goodness-of- 
fit tests that gamma distributions whose parameters can be estimated 
from actual operational data are more precise estimates of the actual 
underlying CONREP service time distributions. Furthermore, it is shown 
that the distribution can be made Erlang by minor adjustments to the 
parameter estimates of the gamma distribution. Such a procedure might 
be desirable for analytical models employing Laplace transforms. 

The data is subjected to linear regression techniques in an 
effort to develop meaningful and accurate functional relationships 
between service time and customer needs. The results indicate that 
the standard error of the estimating relationships would be too large 
to be of any practical use in planning an underway replenishment 
operation. 
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I 



I. INTRODUCTION 



The United States Navy recognizes that one of the most prac- 
tical ways to increase the effectiveness of our naval forces in main- 
taining' superiority at sea is through the application of Underway 
Replenishment (UNREP). The purpose of these operations is to improve 
unit enduranceby supplying needed logistic items to fleet units with 
minimum increase in vulnerability and with minimum diversion from 
performance of combat missions. The resulting ability of ships to 
remain on or near station is equivalent to adding additional combatants. 

The importance of UNREP has been reemphasized within the 
Department of the Navy in that UNREP responsibility is now a 
NAVSHIPSYSCOM designated Project Office (PMS-90) [1], NAVSHIPS 
Instruction 5430.73 of 25 August 1967 is the implementing documeirt. 

Prior to this the responsibility for UNREP was solely a technical one 
and located in various BuShips branches. 

There are basically three ways in which an UNREP may be accom- 
plished. The primary method used is called Connected Replenishment 
(CONREP). This operation consists of intership horizontal transfers 
via rigs connecting the supply ships and the supported units. An 
increasingly important method, which involves the use of UH-34 or 
UH-46 helicopters to deliver certain items, is that of Vertical 
Replenishment (VERTREP). It has been envisioned that VERTREP will 
someday enable all but non-nuclear ships to avoid the inherent vul- 
nerability imposed during CONREP. The third method entails a com- 
bination of CONREP and VERTREP. 

Considerable effort has and is being expended to reduce the 
total time to conduct an UNREP. These efforts include extensive 
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training, improved delivery techniques and equipment, and newer designs 
in supply vessels. Several attempts have been made to model the under- 
way replenishment operations, both from the analytic and computer 
simulation approach. 

McCullough [2] made an analytic approximation of the replenish- 
ment process by using a multi-stage cyclic-queueing model. The model 
considered M supply ships , called stages, placed in series. These 
ships serviced N combatant units, each of which passed by the supply 
ships in succession. In order to achieve a tractable solution, the 
operation was permitted to repeat itself continuously in order to 
study the steady state behavior of the system. A strong assumption 
used in this study was that’ of exponential services times. The most 
basic objection to this approach is that it is not really representa- 
tive of what actually happens during CONREP. First, the cyclic nature 
of this model' is a great simplification of the more complex sequences 
possible in these operations. Secondly, the operation is finite, and 
therefore’ no steady state solution ever exists. The assumption of 
exponential service times in order to make the model more amenable to 
analysis is not valid, a' fact that this paper will attempt to prove. 
Included in McCullough's work, is an incremental ’ time step computer 
simulation- model to test various sequences and to.make comparisons 
with the' results obtained in the analytic model . Here again, a basic 
assumption was the use of specified exponential service times for the 
supply ships, which did not vary between different types of customers. 

Gordon and Copes [3] developed a deterministic model to 
investigate the optimal scheduling-^ of underway replenishment operations 
by considering’ the process as a job scheduling problem. The service 
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times of each combination of supply ship and customer were assumed 
to be deterministic and known, and general expressions were obtained 
for the total time to complete the CONREP and the total waiting time 
of the ships involved. The most serious deficiency of this approach 
is the assumption of fixed service times. It is more reasonable to 
assume that such service times are in fact random variables, and that 
questions relating to the determination of the average time needed to 
complete an underway replenishment operations must take this into 
account. 

Waggoner [4] has developed an analytic model for the case of 
two supply ships and L combatants using queueing theory concepts and 
a random walk in the plane. The distribution of the total time to 
complete the finite operation is obtained in terms of its Laplace 
transform, which although frequently uninvertable, can provide informa- 
tion regarding the mean and variance of CONREP completion times. Here 
again, however, exponential service times were assumed. Waggoner 
concedes the point that although the exponential distribution fre- 
quently fits many realistic queueing models, some doubt exists as to 
its applicability in the case of CONREP operations. 
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II. OBJECTIVES 



The objectives of this paper are twofold: 

1. To determine from operational data more accurate estimates 
of the distribution of service times between various participants in 
the CONREP situation. If this can be accomplished, then more accurate 
results from future analytical and computer simulation models on CONREP 
operations can be expected. 

2. Through the use of scatter diagrams and statistical regres- 
sion techniques, determine if functional relationships exist between 
the amount of a specified type of material received by the combatant 
and the service time necessary to provide this material. 
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IV. GENERAL APPROACH TO CONREP DATA ANALYSIS 



A. DETERMINATION OF CONREP SERVICE TIME DISTRIBUTION AND PARAMETER 
ESTIMATES 

From previous experience, it was felt that service times in a 
CONREP operation would vary in at least two ways. First, the total 
requirement and' replenishment configuration would differ among com- 
batants serviced by the same replenishment ship. Secondly, variation 
would occur among servers replenishing the same type of customers, 
since delivery rates vary from server to server. 

One of the first steps taken to determine the distribution of 
service times wasto'look carefully at the data from several specific 
situations, and from the results find a pattern which might make it 
possible to make a' statement regarding the underlying distribution of 
CONREP service times in' a general sense. 

The three pilot cases chosen were those consisting of replenish- 
ment ship types AOE, AFS, and AE servicing DD (Destroyer) type customers 
These cases were selected primarily due to the large number of sample 
points available, but also because it would permit an evaluation of the 
hypothesis held a priori. 

Frequency histograms were constructed from these sets of data, 
one of which is shown in figure 1. Because of the distinctive skewed 
shapes of these graphs, it was decided to investigate the possibility 
of fitting a "gamma" distribution to each of these situations. The 
gamma distribution is discussed in detail in Appendix A. 

Since the two parameters of the gamma distribution can be 
estimated from the sample mean and variance, these statistics were 
computed (utilizing appropriate computer subroutines). With these 
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SERVICE TIME (MINUTES) 

Figure 1. FREQUENCY HISTOGRAM OF DATA FOR DD (CUSTOMER), AE 



parameters estimated, another subroutine developed for this analysis 
called GAMDIS was utilized to provide the necessary computations neces- 
sary to conduct goodness-of-fit tests with the observed data. Since 
the subroutine GAMDIS is a rather unique program involving numerical 
integration and is not to be found in the usual reference literature, 
a copy is provided in Appendix B. 

The chi-square goodness-of-fit test was used as a measure of 
the appropriateness of fit. Having the expected and observed numbers 
of units possessing service times which fall into "cells" which 
represent an exhaustive and mutually exclusive set of intervals for 
which data is held, the problem is to test the hypothesis 

H: P^-=P^Q (1 = 1 ,2,3,... ,k) 

where p^ is the probability of finding a service time in the i^^ cell, 
p^Q is the observed relative frequency of the actual data falling 
within the same intervals, and k denotes the number of cells. Ten 
minute cells were utilized, with the exception of the last cell, which 
had an interval of time ranging from the end of the final cell having 
a data point within it to plus infinity. The test procedure is to 
calculate the chi-square statistic 

= I (0 -E,)VE, 

1=1 

where 0^. represents the number observed in the i^j^ cell and E^.=np^.^ 
represents the number expected in cell i if H is true, with n equal 
to the sample size. The decision rule is 

reject H If x" i x"(i.„)(k-3)- 

It is necessary to use k-3 degrees of freedom since the two population 
parameter estimates are obtained from the data. 
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The results obtained employing the above procedures on the 

three pilot cases are shown in Table I. In the table, MEAN T denotes 

s 

average service time, VAR. is the sample variance, STD. DEV. is the 
sample standard deviation, r and X are estimates of the parameters of 
the gamma distribution as computed from the data, is the chi-square 
statistic, and is the tabled chi-square value with k degrees 

of freedom, and significance level of a = .05. 



Table I. 

Statistics on Determination of Estimates of Parameters for 
Gamma Distribution for Pilot Cases 



SERVER 


MEAN 

Ts 


VAR. 


STD. 

DEV. 


/N 

r 


X 


X^ 


^\,.95 


COMBINED 


50.2 


462.7 


21.5 


5.4 


.11 


24.5 


35.2 


AOE 


48.0 


308.0 


17.5 


7.5 


.16 


23.3 


26.7 


AE 


50.6 


562.2 


23.7 


4.5 


.09 


21.7 


35.2 


AFS 


51.7 


233.8 


15.0 


11.9 


.23 


8.4 


25.0 



These results indicated that, for the particular cases chosen, suitable 
parameter estimates could be obtained from the data which enabled gamma 
distributions to be fitted to each situation. As indicated, the cal- 
culated value of the chi-square is somewhat less than the tabled value 
in each case. 

As noted in Appendix A, the gamma distribution becomes Erlang 
if the parameter r is a positive integer. Since the Erlang is somewhat 
easier to use in analytic work involving the use of Laplace transforms, 
it was decided to adjust the value of r in each case to the nearest 
integer, and then find the value of X which minimized x^. The results 
of this procedure as applied to the three pilot cases are contained in 
Table II. 
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Table II. 



Statistics on Determination of Estimates of Parameters for 
Gamma Distribution for Pilot Cases, 
f Converted to Nearest Positive Integer 



SERVER 


r 


/s. 

X 


X^ 


y2 

^ .95 


COMBINED 


5.0 


.11 


26.2 


35.2 


AOE 


7.0 


.14 


22.3 


26.7 


AE 


4.0 


.08 


19.4 


35.2 


AFS 


12.0 


.23 


8.3 


25.0 



It is interesting to note that an improvement in was ob- 
tained in all but one case. For this reason, and also due to the 
advantage of working with integer values of f whenever possible, the 
practice was followed throughout the remainder of the effort associated 
with service time distribution computations. Since it was always 
necessary to compute the actual value of the f estimate and its 
associated X^, it was possible to note that the practice improved 
(lowered) the X^ value approximately 60% of the time. It never raised 
the X^ high enough to cause it to fail a goodness-of-fit solely because 
of this adjustment. 

It should be mentioned at this point that no consideration had 
yet been given to the possibility that the distribution of service time 
might also vary depending upon whether the CONREP was conducted during 
daylight or at night. Eventually surmising that this might be the case, 
the data was separated and the analysis rerun. As expected, the average 
service time for night time operations was generally higher, and adjust- 
ments were necessary in the estimated parameters r and a in each case. 
Therefore all remaining calculations were conducted on data which had 
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been carefully separated between day and night operations. Fortunately, 
all data sources clearly indicate this important factor. 

B. LINEAR REGRESSION MODEL TO DETERMINE THE RELATIONSHIP BETWEEN CONREP 
SERVICE TIME AND AMOUNT OF MATERIAL TRANSFERRED 

Prior to undertaking this part of the data analysis, it was felt 
that nearly 1 inear relationships existed between customer needs and the 
times necessary to provide them. If* these relationships could be found, 
they could be a valuable aid to those responsible for planning an UNREP. 
For example, if this technique could be utilized to determine the total 
time of an UNREP, then this knowledge would aid planners and partici- 
pants alike since more effective use of combatants and Underway Replen- 
ishment Groups (URGs) would result. 

In order to provide a visual aid in determining what relation- 
ships actually existed between service times and requirements, scatter 
diagrams were produced using BIMED 02D [10], a statistical subroutine 
contained' in the function library of the IBM 360 computer facility at 
the Naval Postgraduate School. The cases chosen as representative of 
the many possible combinations were those involving AOs and AEs as 
replenishment vessels and combatant types CVA, CLG, and DD. This 
selection was chosen since AOs provide an assortment of liquid fuels, 
while AEs supply only dry material, namely, ammunition and explosives. 

Before discussing the scatter diagrams further, it is neces- 
sary to digress for the moment to describe the unique problems 
presented by the AO server situation. These supply vessels have the 
capability of supplying three different fuels, namely. Navy Standard 
Fuel Oil (NSFO), Diesel Fuel (JP-5), and Aviation Gasoline (AVGAS). 

One customer, the conventionally powered aircraft carrier, frequently 
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takes all three types of fuel simultaneously through separate hoses. 
Other customers, say helicopter equipped, conventionally powered 
cruisers and destroyers, will take NSFO or JP-5 and perhaps AVGAS. 

As can be seen on the' scatter diagrams, there is a wide range of fuel 
type and requirements' combinations. Thus the problem of finding a 
single, controlling independent variable for the AO server situation 
would' appear to have no solution. The alternative is to attempt to 
find a method of predicting the- response based on the magnitudes of 
a number of variables. 

The scatter diagrams for 16 situations involving the afore- 
mentioned participants are shown in figures 3 through 18 in Appendix C. 

In the special case of the AO serving CVA/CVS ships, each of 
the independent variables, barrels of NSFO, JP-5, and AVGAS, were 
plotted individually against total CONREP service time. In all other 
cases, only the dominant independent variable was plotted with time. 

All daylight CONREP situations are plotted separately from night 
operations. 

A discussion of ' the scatter diagrams is contained in the 
next section. 

The technique utilized to discover the functional relation- 
ships between needs and service times was that of linear regression. 

In the AO-CVA/CVS cases, usually involving three independent variables, 
multiple linear regression was applied. The aim was to determine if 
a regression model of the form 

T=a+6X +3X +GX +E 

S 1 1 2 2 3 3 

could be calculated and applied in this case, where T^ represents the 
response dependent variable, total service time, and X^, X^, and X^ 
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represent the independent variables, barrels of NSFO, JP-5, and AVGAS 
respectively. E denotes the error term in the model. 

Since the remainder of the cases involve only one independent 
variable, the attempt was to find a simple linear regression model of 
the form 

= a + 3X + E 

where X represents barrels of NSFO for the other AO serving situations, 
and short tons of explosive material for all AE serving situations. 

BIMED 03R, a multiple regression subroutine contained in the 
function library of the IBM 360 computer facility at the Naval Post- 
graduate School, was utilized for all regression work. 

No attempt was made to fit higher order polynomials to the 
CONREP data. 
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V. RESULTS OF THE DATA ANALYSIS 



A. DETERMINATION OF CONREP SERVICE TIME DISTRIBUTIONS AND PARAMETER 
ESTIMATES 

The techniques described in Section IV for the determination 
of the parameters of the gamma distribution from the data were applied 
to a total of 29 out of a possible 36 separate cases, consisting of the 
combinations of' AO, AE, AOE, AFS, AF, and AKS serving CVA/CVS, CLG/DLG/ 
DDG, and DD type customers, for day and night operations separately. 

The results of this determination are contained in Tables III 
through VIII. The following abbreviations are used in the tables. 

N = total number of sample points considered or 
available. A maximum of 30 was utilized to 
reduce size of problem. 

T^ = average CONREP service time of sample. 

Ofs = sample standard deviation of CONREP service time. 

r = estimate of parameter r for gamma distribution 
to nearest integer. 

X = estimate of parameter X for gamma distribution. 

= chi-square statistic for goodness of fit. 

k = degrees of freedom. 

X^. Q 5 = tabled value of chi-square statistic with k degrees 

“ of freedom and significant level a = .05. 

Statistics were not computed for situations in which eight or 
less data points were held. 

Only two cases, AO-DD (Night) and AFS-CVA/CVS (Day) failed the 
chi-square goodness-of-fit tests at the a = .05 level. Both would have 
passed at the .03 level. 

Eight out of 29 cases had estimates of the r parameter of 9 or 
greater, indicating that these cases could be estimated by a normal 
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population with parameters of mean and standard deviation indicated 
in the tables. 





Statistics 


Table III. 
on AO Serving 


Situation 








CVA/CVS 


CLG/DLG/DDG 


DD 






DAY 


NIGHT 


DAY 


NIGHT 


DAY 


NIGHT 


N 


30 


22 


30 


30 


30 


30 


Ts 


114.3 


123.4 


58.0 


72.1 


81.2 


78.3 


/s 

^ts 


34.1 


30.4 


21.8 


22.1 


41.1 


29.6 


/V. 

r 


11.0 


16.0 


7.0 


11.0 


4.0 


7.0 


✓s 

X 


.09 


.13 


.11 


.15 


.05 


.09 


X2 


20.5 


15.9 


20.2 


9.9 


27.3 


28.4* 


k 


19 


17 


13 


12 


17 


17 


y2 

^ k,.95 


30.1 


27.6 


22.4 


21.0 


27.6 


27.6 




Statistics 


Table IV. 
on AE Serving 


Situation 








CVA/CVS 


CLG/DLG/DDG 


DD 






DAY 


NIGHT 


DAY 


NIGHT 


DAY 


NIGHT 


N 


30 


30 


30 


30 


30 


30 


Ts 


98.6 


161.2 


71.8 


181.3 


41.2 


83.3 


^ts 


61.2 


52.8 


51.5 


123.0 


23.7 


36.1 


r 


3.0 


9.0 


2.0 


2.0 


3.0 


5.0 


/s 

X 


.03 


.06 


.02 


.01 


.06 


.07 


X2 


25.8 


27.6 


34.7 


62.4 


12.0 


25.3 


k 


26 


27 


26 


51 


10 


16 


>=\,.95 


38.9 


40.1 


38.9 


68.7 


18.3 


26.3 
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Table V. 



Statistics on AOE Serving Situation 





CVA/CVS 


CLG/DLG/DDG 


DD 






DAY 


NIGHT 


DAY 


NIGHT 


DAY 


NIGHT 


N 


19 


5 


23 


3 


30 


11 


Ts 


141.7 


1 

1 

1 

t 


81.2 


t 

1 

1 

t 


51.6 


70.9 




58.9 


1 

<c 
1 — 


41.1 


1 

c 

t— 


18.3 


35.7 


ts 

y\ 

r 


6.0 


c 

o 


4.0 


<x, 

Q 

I— 


8.0 


4.0 


X 


.03 


'Zl 

LU 

1—4 

o 


.05 


Z1 

LU 

H- « 

O 


.17 


.05 




28.0 


1—4 

U. 

u. 


27.1 


H- 1 

U. 

LU 

=D 


12.7 


18.8 


k 


31 


oo 

Z1 

1 

1 

1 


17 


CO 

•— 1 
1 
1 


8 


13 


^ k,.95 


45.0 


1 

1 

1 


27.6 


1 

t 

i 


15.5 


22.4 
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on 


Table VI. 

AFS Serving Situation 






CVA/CVS 




CLG/DLG/DDG 


DD 






DAY NIGHT DAY 


NIGHT 


DAY 


NIGHT 


N 


13 


3 


14 


13 


30 


24 


Ts 


145.8 


1 

1 

1 

1 


108.3 


71.5 


52.4 


58.5 


/\ 

<^ts 


61.6 


1 

<c 

§ 


60.0 


30.3 


12.2 


21.3 


/V 

r 


5.0 


LU 


3.0 


5.0 


18.0 


7.0 


✓N 

X 


.03 


fc— c 

O 

1—4 

LU 


.02 


.07 


.35 


.12 




46.5* 


Lu 

ID 

LO 


27.3 


22.0 


6.8 


4.3 


k 


30 


HH 

1 

1 


26 


12 


6 


9 




43.8 


1 

1 

1 


38.9 


21.0 


12.6 


16.9 
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Table VII. 

Statistics on AF Serving Situation 





CVA/CVS 




CLG/DLG/DDG 


DD 






DAY 


NIGHT 


DAY 


NIGHT 


DAY 


NIGHT 


N 


9 


6 


16 


9 


30 


20 




113.4 


1 

1 

1 

1 


61.9 


69.2 


55.2 


54.3 


/s 

^ts 


45.3 


1 

< 

h- 

c 

Q 


29.2 


27.7 


17.4 


21.1 


r 


6.0 


1 — 
LU 


4.0 


6.0 


10.0 


6.0 


X 


.05 


1— 1 
O 
»— 1 
Lu 


.06 


.08 


.18 


.11 


X2 


25.9 


Ll. 

ZD 

CO 


16.7 


9.1 


8.2 


11.1 


k 


17 


1 — 1 

1 

1 


13 


11 


7 


9 
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Statistics 


Table VIII. 
on AKS Serving 


Situation 






CVA/CVS 




CLG/DLG/DDG 




DD 




DAY 


NIGHT 


DAY 


NIGHT 


DAY 


NIGHT 


N 


3 


6 


20 


5 


30 


19 


Ts 


1 

1 


1 

1 

1 


41.3 


1 

1 

1 

i 


36.2 


45.9 


/s 

^ts 


1 

1 

1 

< 

h- 


1 

1 

< 
1 — 
< 


13.5 


1 

< 

1— 

< 

o 


11.9 


20.4 


/s 

r 


< 

Q 

h— 


Q 

1— 


9.0 


I— 

LU 


9.0 


5.0 


/N 

X 


LU 

t— 1 

O 

1 — i 


LU 

1 — i 

O 
1— 1 
Li. 


.21 


1 — i 

O 

1 — i 

Li. 

Li. 


.24 


.10 




Li. 

Li. 

ZD 

CO 


Li. 

ZD 

CO 


2.7 


ZD 

CO 

1— 1 


6.3 


8.0 


k 


1— 1 
1 


t—t 

1 

1 


6 


1 

1 

1 


6 


9 


’^^,.95 


1 

1 

1 

1 


1 

1 

1 


12.6 


1 

1 


12.6 


16.9 
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B. LINEAR REGRESSION MODEL TO DETERMINE THE RELATIONSHIP BETWEEN 
CONREP SERVICE TIME AND AMOUNT OF MATERIAL TRANSFERRED . 

With few exceptions the scatter diagrams in Appendix C tend to 
discourage the notion that meaningful, accurate linear relationships 
exist between combatant ship requirements and service time. This is 
especially true in the case of the AO servicing various units. There 
is a definite improvement in those situations involving the AE as the 
server. On several of these diagrams, it is possible to fair in a 
line which would serve as a linear model. It is questionable from 
the results' of the scatter diagrams that the variability around a 
single' regression' line could be low enough to be acceptable. In order 
to determine this vari ability' more precisely, the multiple regression 
techniques mentioned in Section IV were employed. The results of these 
efforts are shown in Table IX, for the AO and AE server situation only. 

As an example, this table indicates that a linear regression 
model for predicting the amount of time necessary for an AO to supply 
quantities of NSFO, JP-5, and AVGAS could be computed using the 
following formula: 

T, = 68.1 + .0050 X + .0014 X + .0062 X (minutes) 

s 1 2 3 

where X denotes barrels of NSFO, X barrels of JP-5, and X barrels 

1 2 3 

of AVGAS required. The standard error (i.e., the standard deviation) 
of the estimate is 27.0 minutes. Three cases were selected from the 
dependent sample to illustrate the accuracy of this procedure. The 
results were as follows: 



NSFO 

(BBLS) 


JP-5 

(BBLS) 


AVGAS 

(BBLS) 


Ts 

CALCULATED 


Ts 

ACTUAL 


DIFFERENCE 

OBSERVED 


5681 


3535 


1212 


108 


96 


+12 


4139 


12452 


0 


106 


99 


+ 7 


4595 


0 


0 


90 


158 


-68 
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If this sample can be considered representative, it illustrates quite 
clearly that magnitudes of the errors that might be expected. 

As expected and in agreement with the scatter diagram, the 
standard estimate errors for the cases involving the AE as server were 
less than those involving AOs in four out of six cases. 
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Table IX, 

Results of Linear Regression Analysis on AO and AE 
Server Situations 



AO-CVA/CVS (Day) 

Sample size 

Percent variance explained 
Multiple correlation coefficient 
Standard error of the estimate 
a - intercept value 

Mean Std. 

Dependent variables (BBLS) Dev. 


45 

36.4 

.604 

27.0 min. 

68.1 min. 

Reg. Std. Err. Partial Corr. 
Coe. Reg. Coe. Coefficients 


NSFO 


6577.9 


3460.5 


0.0050 0.00125 


0.53 


JP-5 


6267.8 


4638.2 


0.0014 0.00093 


0.22 


AVGAS 


273.2 


655.3 


0.0062 0.00639 


0.15 


AO-CVA/CVS (Night) 












Sample size 






20 






Percent variance explained 




22.6 






Multiple correlation 


coefficient 


.4757 






Standard error of the estimate 


32.8 min. 






a - intercept value 






54.9 min. 








Mean 


Std. 


Reg. 


Std. £irr. 


Partial Corr. 


Dependent variables 


(BBLS) 


Dev. 


Coe. 


Reg. Coe. 


Coefficients 


NSFO 


6939.5 


4579.7 


0.0027 


0.00173 


0.36 


JP-5 


5454.8 


3136.9 


-0.0013 


0.00251 


-0.13 


AVGAS 


258.3 


514.4 


0.0212 


0.01489 


0.34 


AO-CLG/DLG/DDG (Day) 












Sample size 






51 






Percent variance explained 




43.2 






Standard error of the estimate 


14.4 min. 






a - intercept value 






24.8 min. 








Mean 


Std. 


Reg. 


Std. Err. 


Partial Corr. 


Dependent variable 


(BBLS) 


Dev. 


Coe. 


Reg. Coe. 


Coefficient 


NSFO 


2025.1 


879.0 


0.0141 


0.00230 


0.66 


AO-CLG/DLG/DDG (Night) 










Sample size 






49 






Percent variance explained 




10.0 






Standard error of the estimate 


24.0 min. 






a - intercept value 






56.5 min. 








Mean 


Std. 


Reg. 


Std. Err. 


Partial Corr. 


Dependent variable 


(BBLS) 


Dev. 


Coe. 


Reg. Coe. 


Coefficient 


NSFO 


2111.0 


1470.9 


0.0054 


0.00236 


0.32 
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Table IX. (cont.) 



AO- DP (Day) 



Sample size 99 

Percent variance explained 22.0 

Standard error of the estimate 30.8 min. 

a - intercept value 25.3 min. 

Mean Std. Reg. 

Dependent variable (BBLS) Dev. Coe. 



NSFO 1320.9 533.2 0.0305 



Std. Err. 
Reg. Coe. 

0.00583 



AO-DD (Night) 

Sample size 

Percent variance explained 
Standard error of the estimate 
a - intercept value 



40 

0.02 

60.8 min. 

80.9 min. 



Dependent variable 
NSFO 



Mean 

(BBLS) 

1490.7 



Std. Reg. Std. Err. 

Dev. Coe. Reg. Coe. 

449.0 0.0020 0.02164 



AE-CVA/CVS (Day) 

Sample size 

Percent variance explained 
Standard error of the estimate 
a - intercept value 



30 

85.5 

23.4 min. 
-30.4 min. 



Mean ^ Std. 
Dependent variable (S.T. ) Dev. 



Reg. Std. Err. 

Coe. Reg. Coe. 



AE Material 



207.2 126.9 0.4332 0.03419 



AE-CVA/CVS (Night) 

Sample size 

Percent variance explained 
Standard error of the estimate 
a - intercept value 



30 

68.2 

31.9 min. 
48.6 min. 



Mean Std. 
Dependent variable (S.T. ) Dev. 



Reg. Std. Err. 

Coe. Reg. Coe. 



AE Material 



285.6 122.0 0.376 0.04861 



* S.T. denotes Short Tons. 



Partial Corr. 
Coefficient 

0.47 



Partial Corr. 
Coefficient 

0.0151 



Partial Corr. 
Coefficient 

0.923 



Partial Corr. 
Coefficient 

0.826 
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Table IX. (cont.) 



AE-CLG/DLG/DD6 (Day) 

Sample size 

Percent variance explained 
Standard error of the estimate 
a - intercept value 



Dependent’ variable 


Mean 

(S.T.) 


Std. 

Dev. 


AE Material 


50.7 


53.8 


AE-CLG/DLG/DDG (Night) 




Sample size 

Percent variance explained 
Standard error of the estimate 
a - intercept value 


Dependent variable 


Mean 

(S.T.) 


Std. 

Dev. 


AE Material 


90.8 


85.7 


AE-DD (Day) 






Sample size 

Percent variance explained 
Standard error of the estimate 
a - intercept value 


Dependent variable 


Mean 

(S.T.) 


Std. 

Dev. 


AE Material 


18.0 


16.0 


AE-DD (Night) 






Sample size 

Percent variance explained 
Standard error of the estimate 
a - intercept value 


Dependent variable 


Mean 

(S.T.) 


Std. 

Dev. 


AE Material 


27.9 


21.4 



.7 

.7 min. 
.3 min. 



Reg. 

Coe. 


Std. Err. 
Reg. Coe. 


Partial Corr. 
Coefficient 


0.965 


0.07144 


0.93 


.8 

.5 min. 
.4 min. 






Reg. 

Coe. 


Std. Err. 
Reg. Coe. 


Partial Corr. 
Coefficient 


1.337 


• C.1489 


.882 


.5 

.3 min. 
.7 min. 






Reg. 
. Coe. 


Std. Err. 
Reg. Coe. 


Partial Corr. 
Coefficient 


0.568 


0.0543 


0.892 


.3 

.1 min. 
.4 min. 






Reg. 

Coe. 


Std. Err. 
Reg. Coe. 


Partial Corr. 
Coefficient 


1.67 


0.34 


0.68 



30 

86 

20 

80 

25 

77 

62 

67 

30 

79 

7 

23 

30 

46 

39 

36 



32 



VI. CONCLUSIONS 



It has been shown that the gamma distribution, with parameters 
determined from operational data, can be used to describe the probability 
distribution of CONREP service times for a large number of practical 
situations. With minor adjustments in the gamma parameter estimates, 
the Erlang distribution, with its increased flexibility in analytical 
models employing Laplace transforms, can be used effectively for the 
same purpose. 

It is conceivable that previous work in analytical and computer 
simulation studies which assumed exponentially distributed CONREP service 
times can be extended to incorporate this new information, with a result- 
ant increase in the accuracy of these models in depicting actual under- 
way replenishment operations. 

The second part of the analysis indicates that, with few excep- 
tions, estimates of service time utilizing linear regression techniques 
would be too inaccurate, as indicated by the standard error of the 
estimates derived from the data, to serve the purposes for which such 
a model might be used. Although certain cases appeared to be quite 
accurate, the overall results discourage enthusiasm for further effort 
along these lines. It is noted that some cases appear to lend them- 
selves to non-linear estimates, such as might be obtained through poly- 
nomial regression, since some curvature can be detected in several of 
the scatter diagrams. A larger data base than that which was available 
for this study would be necessary to make meaningful conclusions along 
these lines. 
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APPENDIX A 

The Gamma Distribution 



If a process consists of r successive events and if the total 
elapsed time of this process can be regarded as the sum of r independ- 
ent exponential variates each with parameter X, the probability 
distribution of this sum will be a gamma distribution with parameters 
X and r. A form of the gamma distribution may be fitted to many 
positively skewed distributions of statistical data. 

The gamma distribution is described by the following density 



function 

(r-1) -XX 

g(x;r,x) (p_i) j 

where X > 0, r > 0, and x are non-negative. 

The sum of r (where r is a positive integer) exponential 
variables each having the same parameter X is also called an Erlang 
distribution. Mathematically, the Erlang distribution is a convolu- 
tion of r exponential variates, i.e., the distribution of the sum of 
r exponential variables. Because of the restriction on r, a process 
which can be represented by an Erlang distribution is easier to work 
with in an analytic sense, especially if Laplace transform techniques 
are utilized. 

The expected value and variance of the gamma distribution are 



given by 
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From these relationships it follows that 



r = 



. (EX)= 



vx 



X = ^ 



EX 

VX 



Therefore, a gamma distribution, like a normal distribution, 
can be specified by its mean and its variance. 

If r = 1, the gamma distribution is identical to the exponen- 
tial distribution 



f(x) = Xe"^^ X ^ 0. 

As r increases, the gamma distribution approaches a normal 
distribution asymptotically, r is often referred to as the "shape" 
parameter, while X is known as the "scale" parameter. Figure 2 
illustrates the graphs of g(x;r,X) for several values of r and X = 1. 

In the gamma distribution the ratio of the standard deviation (the 
square root of the variance) to the mean (expected value) is the 
reciprocal of the square root of r. When this ratio is 1/3 or 
smaller (i.e., when r is 9 or greater), the gamma distribution closely 
approximates a normal distribution with the same mean and standard 
deviation. 

The cumulative distribution function does not exist in explicit 
form for the gamma distribution, although the values of the so-called 
incomplete gamma function have been tabled by Pearson [7]. The computer 
subroutine entitled GAMDIS was constructed for and used in this analy- 
sis. The subroutine employs numerical integration techniques involving 
the trapezoidal rule [8]. GAMDIS has the function of computing proba- 
bilities and expected numbers necessary in making tests of hypotheses 
in chi-square goodness-of-fit tests. Because of the potential usefulness 
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of this subroutine to others engaged in future efforts of this nature, 
a copy of the subroutine, written in FORTRAN with explanatory comment 
cards, is included as Appendix B. 
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Figure 2. GAMMA PROBABILITY DISTRIBUTION IN THE SPECIAL CASE WHERE X 



APPENDIX B 
Subroutine GAMDIS 



Subroutine GAMDIS is called by specifying the estimates of the 
gamma distribution, R and XLAMDA. These estimates are computed from 
the sample mean and sample variance of the service times under con- 
sideration in accordance with the following relationships: 

R = p = (sample mean)^ 
sample variance 

XLAMDA = X = sample mean 
XLAMDA A variance 

Other calling parameters for the subroutine include NBU, the number of 
units in the sample, and NBCL, the number of cells for which expected 
numbers of units having appropriate service times are to be calculated 
and labeled by vector EXPECT. 

The subroutine performs the following computations: 

1. Defines the limit of the total area for which numerical 
integration will be utilized (LIMIT = NBCL*10). 

2 . Defines constants based on given values of R and XLAMDA 
which will be used in specifying the gamma density function (F) used 
in numerical integration. 

3. Sets up subintervals necessary to obtain an accurate 
measurement of the probabilities associated with each interval. 

4. Compute the probabilities mentioned in 3., (SUM2), plus 
the remainder of probability associated with the last interval (SUM4). 

5. Computes the expected number of units for each cell using 



NBU, SUM2, and SUM4. 

6. Returns the computations of 5. to the main program via 
vector EXPECT. 
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StlRROKTINE GAMOr, ( R , XL A ^'^)A , N BU , NRf. L , FX F EC T,» 

IT T yr. n. v. i^ ^ 

r I MENSTON A (5i : ) ,ERAP ( '^CC ) » AJ( ICOl ) ,FXCECT( 5'' > ' 

^ P=SHAPE PAR4MFTFP OF GAMMA OISTPIBUTION 

C XLAMDA = scale Pi<?AMETFR 

C MRU = total number OF SAMPLE UNITS 

C NRCL = TOTAL NUMBER OF 'CELLS* TO RE CONRIDFPFO 

C EXPECT = VECTOR OF EXPECTED NUMBER UNITS IN EACH 'CELL' Rf'^‘= 

C CN GIVEN R ANO XLA'^OA 

LIMIT = NRCL’*1'“ 

PLAMBA=XLAMnA tP 

N=10 

Nl = l^ 

AN = N 

AN1=N1 

GX=GAMMA(R) 

CONS=PL AMBA/GX 
L = 0 

SUMR=C 

ro 15 jj = n, LIMIT, u 
^^0 = 0 0 
PLIM=PJ 
/ LIM = P 

XX = ( BLI "^-AL /AN 

A(1)=ALIM 

NK=N+1 

no ?r K=2,NK 
A <=K-1 

"*0 A ( K) =AL IM + XX* AK 
5UM2=r 
ro 5^ K=1 ,N 
SUM1=0 
NK=N1+1 
no Af 1=?. ,NK 
AI=I-I 
XXX=XX/AN1 
AJ( I )=A(K )+XXX» AI 
IF (AJ( I ).LE.->,C » goth 4C 
7 - YV = AJ{ I ) 

F = EXP(-XLAMr)A>-YY) 

C F = VALUE OF GAMMA DENSITY FUNCTION AT X=YY 

F = CnNS»YY = -<‘ ( P-1. ) 'F 
SUM1 = SUM14-F 
4'. CONTINUE 

C FRAR(K) = AVERAGE VALUE OF F IN SUBINTFRVAL 

FRAR(K)=SUM1/(AN1+1. ) 

C T1 = AREA CF SUBINTFRVAL 

■*■)=( A ( K+l )-M K ) )*FBAP ( K ) 

C SUM2 = probability ASSOCIATED WITH THIS 'CELL* 

'^UM2 = SUM2+T1 
5^ CON-'-INUE 
L = L + 1 

FXPECT(L)=SUN2»NBU 

C SUM3 = CUMULATIVE PROBABILITY THROUGH LA^T COMPUTATION 

SUM3 = SUM3*SUM2 ^ 

15 CONTINUE 

C 5UM4 = PROBABILITY ''CSOCIATEO WITH LAST CFLL 

SIIM4 = 1,C-SUM’ 

EXPECT! L + 1 ) =SUM4- NBII 

RETURN 

END 
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APPENDIX C 

Scatter Diagrams of Operational CONREP Data for Situations 
Involving AO, AE Servers, CVA/CVS, CLG, and DD Customers 
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FIGURE 3. BI^'EC C2C SCATTER CIAGBAF* FOR AC-CVA/CVS (0AY» SITUATION, NFSO VS SERVICE TIME 
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FIGURE 4. 8IMEC C2D SCATTER CIAGRAH FOR AC-CVA/CVS (CAYJ SITUATICNf JP-5 VS SERVICE TIME 
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FIGURE 5. BIMEC 02T SCATTER DIAGRAM FOR AC-CVA/CVS (DAY) SITUATICN* AVGAS VS SERVICE TIME 
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FIGURE 6. BIMEC C2D SCATTER DIAGRAM FOR AC-CVA/CVS (MGHTI SITUATION, NSFO VS SERVICE TIME 
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FIGURE 7. BP'EC C2C SCATTER DIAGRAM fCR AC-CVA/CVS IMGHTI SITUATICN, JP-5 VS SERVICE TIME 
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FIGURES. BIHEC 02 C SCATTER DIAGRAM FOR AC-CVA/CVS (NIGHII SITUATION, AVGAS VS SERVICE TIMp 
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FIGURE 9. BII*ED 020 SCATTER DIAGRAM FOR AO-CIG/OLG (DAY! SITUATION, NSFO VS SERVICE TI*"E 
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FIGURE 10. BIHEC 02D SCATTER ClAGRAH FCR AC-CLG/DLG IMGHTI SITUATICNt NSFO VS SERVICE TIME 
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FIGURE 11. AIMED 02D SCATTER CIAGRA- FOR AC-CC ICAY) SITUATION, NFSC VS SERVICE TIME 
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FIGURE 12. BII'EC 02C SCATTER DIAGRAM FCP AC-CC (MGHTI SITUATION 



VARIABLE SHORT TENS OF AE MATERIAL 
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FIGURE 13. BIHEO 02C SCATTER CIAGRAM FOR AE-CVA/CVS (DAY> SITUATICN, AE MATERIAL VS SERVICE TIME 
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FIGURE 14. BIMEC 020 SCATTER CIAGRAM FOR AE-CVA/CVS (NIGHT! SITUATICN, AE MATERIAL VS SERVICE TIME 
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FIGURE 15. BIHEO 02C SCATTER ClAGRAM pOR AE-CLG/CLG (CAYI SITUATION, AE MATERIAL VS SERVICE TIME 
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FIGURE 17. 
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FIGURE 18. BIMEC C20 SCATTER DIAGRAM FOP AE-CC (MGHTI SITUATION, AE MATERIAL VS SERVICE Time 
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